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Abstract. The formation of planetesimals in protoplanetary disks due to collisional sticking of smaller dust aggre- 
gates has to face at least two severe obstacles, namely the rapid loss of material due to radial inward drift and 
particle fragmentation due to destructive collisions. We present a scenario to circumvent these two hurdles. Our 
dust evolution model involves two main mechanisms. First, we consider a disk with a dead zone. In an almost 
laminar region close to the midplane, the relative velocities of the turbulent particles are comparatively small, 
which decreases the probability of destructive particle collisions. Second, turbulence is not the only source of vi- 
olent relative particle velocities, because high radial drift speeds can also lead to boulder fragmentation. For this 
reason, we focus additionally on the snow line. Evaporation fronts can be associated with gas pressure maxima 
in which radial drift basically vanishes. This implies that particle fragmentation becomes even less likely. Our 
simulation results suggest that particles can overcome the fragmentation barrier. We find that boulders of several 
10 2 m can form within only a few thousand years. 
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1. Introduction 

The process of planetary formation is widely accepted 
to be initiated in circumstellar disks by the colli- 
sional coagula t ion o f dust grains to larger objects 
(jHenning et all l2006t iNatta et alj l2007ri . Evidence of 



insta b ility dGoldreich fc Ward! Il973t IWeidenschillind 
19791: lYoudin fc Shul 12061" ISchrapler fc Hennind 120041 ; 
Garaud fc LinT hoO^ - an issue that is still highly debated. 



Due to the overwhelming amount of observational data 
for protoplanetary disks available, models o f dust particle 
grow t h have attracted more at t ention dTanaka et al ' 



2005[ iDullemond fc Dominikl 120051: iNomura fc Nakagawa 



grain growth to particle sizes beyond those typically 
observed in the interstellar medium is provided by 



2006; iD'Alessio et all 120061; lOrmel et al.1 l2007t iCiesla 
20071 ). All of these particle evolution models can ex- 



plain observational disk features su ch as, for example 



(Bouwman et al.ll200lt 


van Boekel et al. 2003; Apai et al. 


2004 Kessler-Silacci et al. 


2007: 


Sicilia-Aguilar et al. 


20071) . Analyses of millimeter interferometrv data of disks 



the group I/II disk classif ication (jMeeus et al.1 12001 



Dullemond fc Dominikl l2004f ) or the disappearance of 



infrared excess in t he spectra of disks with ages higher 



2001; iCarpenter et al 



grains exist with 



implies that large populations of dust 



particle radii up to several millimeters ( Testi et al. 2003; 
Wilner et al]|2005t iRodmann et al.ll2006l ). 

Models of dust particle growth in protoplanetary 
disks predict coagulation, which supports the idea of 
planetesimal formation by particle hit and stick mech- 
anisms. Numerical calculations support the formation 
of cm-sized particles within a few 10 3 orbital timescale s 
in the inner parts of the disk ( WeidenschiHind 1980f ). 
After particles have grown to larger sizes, they can also 
be affected by other processes described by alternative 
planet formation scenarios. For example, large grains 
experience vertical settling, which invol ves the forma- 
tion of a dense midplane dust layer (jGaraud fc Linl 
200J). This layer can be affected by gravitational 



than a few Myrs (jHaisch et al 
120051) . although, there are still many open observational 
riddles that these models will hopefully answer in the 
near future. 

Models of protoplanetary disks do not only help 
to interpret observational data; they also iden tify the 
serious obstacles to plan etesimal formation (jYoudin 



2004 iDominik et ail l2007h . One of these obstacles is 

Due to the high particle rel- 
that can reach 100 m/s 



particle fragmentation. 

ative velocities in disks? 

(IWeidenschillind Il977t lOrmel fc Ci 



20071 ). particle 



collisions lead to particle destruct i on instead of par- 



ticle growth dBlum fc Wurml |2000l : iBrauer et all 12008 



Johansen et al. 20081 ). Depending on the disk model, dust 



particle growth is inhibited significantly by high-speed im- 
pacts of particles around a meter size. Even for times as 



2 
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long as 1 Myr, solid particles are unable to overcome the 
fragmentation barrier and the formation of planetesimals 
as precursors for Earth-like or Jovian planets, therefore, 
poses a major problem. Despite this theoretically pre- 
dicted inability for planetesimals to form, the orbital decay 
of dust particles can, nevertheless, reproduce the observed 
evolu t ion in the spectral ener gy distribution (Hai sch et al.l 



200 It ICarpenter et al.l 120051 ) . The presence of a growth 



barrier is, therefore, not in contradiction with observa- 
tions. 

Several mechanisms in prot oplanetary dis k s can 



produce planetesimal formation. iJohansen et al.l (|2007f ) 



showed that the non-linear feedback of dust onto gas 
can lead to the rapid formation of gravitationally bound 
chimps of dust, which subsequently form Ceres-size bod- 
ies. The dust particles, however, must have already 
grown to some meters in size before this scenario can 
take place. Another possibility for solving the forma- 
tion problem is particle trapping in gas pressure maxima 
(jBarge fc Sommerialll995tlKlahr fc Henning|ll997l ). In gas 
pressure bumps, relative particle velocities are substan- 
tially lower and, hence, the collision between large bodies 
is more likely to lead to particle growth than to particle 
disruption. 

Interesting means of gas press ure maxima and dus t 
particle trapping was presented by iKretke fc Linl (|2007l ) . 
This particle retention mechanism requires the presence 
of an evaporation front, for example the snow line, which 
acts in the following way. As we travel through the snow 
line in a direction away from the central star, the dust- 
to-gas ratio increases suddenly. This jump in dust density 
affects the strength of the magneto-rotational turbulence, 
since the amount of free electrons in the disks strongly de- 
pends on the dust density. With increasing d ust density. 



the amount of turbulence in the disk decreases ( Sano et al 



2000; Ilgnci fc Nelsod f2006). We assume a constant mass 
accretion rate throughout the disk. The gas surface den- 
sity then has to be higher in low than in high disk tur- 
bulent regions for the gas mass accretion to remain at 
a similar value. The sudden increase in the dust density 
could therefore also reproduce a su dden increase in gas 
densities. For certain accretion rates, IKretke fe Lml (|2007l ) 
found the occurrence of a local gas density bump in which 
solid particles tend to accumulate. In their simulations, 
dust particle retention reproduced very high surface dust 
densities of the order of several 10 3 g/cm 2 , which again 
raises concern about gravitational instability. 

However, even for a particle retention mechanism sim- 
ilar to that described above, the growth of solid mate- 
rial towards larger sizes remains an open issue. Gas pres- 
sure maxima may decrease relative radial particle veloci- 
ties triggering coagulation; the most severe reason for vi- 
olent particle fragmentation is not however radial drift 
but the turbulent nature of the protostellar disk itself. 
Relative turbulent particle velocities are appr oximately 
10 m/s (|Ormel fc Cuzzill2007t IVolk et al.lll980l) and par- 
ticle sticking at the s e high speeds is extremely unlikely. 
Nevertheless, Cieslal ( 2007 ) found that particles can grow 



to several 10 m in radius if we consider layered MM active 
disks. Under certain circumstances, MRI is only active in 
the upper layers of the disk, while the disk midplane is 
almost laminar. Since most of the larger grains are lo- 
cated around the midplane where the disk is quiescent, 
dust growth is not inhibited by the high speed collisions 
that produce particle fragmentation. 

We combine three ingredients for a planetesimal for- 
mation model. We consider dust par ticle growth, parti- 
cle fragmentation and radial motion (IBrauer et "al. 20081 ) 
around the snow line (jKretke fc Linl 12 007) in a layered 
MRI active protoplanetary disk ( Cieslall2007j ). The inclu- 



sion of a snow line in our simulations provides a parti- 
cle retention mechanism which almost erases the radial 
drift velocities favouring particle growth. We adopt a lay- 
ered MRI-drivcn disk to circumvent particle fragmenta- 
tion in the midplane of the disk due to turbulent motions 
of the gas. We investigate whether solid particles can over- 
come the fragmentation barrier and produce larger ob- 
jects, which represent the possible precursors to the plan- 
ets. We study the influence of two parameters, namely the 
gas accretion rate M acc and the critical threshold velocity 
for fragmentation Vf . 

2. Model 

We assume a background gas disk which is in a steady 
state and, hence, does not change with time. The calcu- 
lation of the gas d ensities around the snow line are de- 
scribed in detail in IKretke fc Linl ( 2007 ). We adopt all of 
the parameter values presented in that paper apart from 
the residual turbulence value around the midplane ao, 
which we set to be 10~ 5 corresp o nding to self- induced tur- 
bulcncc (Wcidcnschilling 119791 : IWeidenschilling fc Cuzzil 
19931 ). The water evaporation front, which is the impor- 



tant element in our model, is located at 3 AU (|Lecar et al 
20061 ) . The gas surface densities that we adopt in our sim- 



ulations are shown in Fig. [T] for different accretion rates 

Macc 

We focus on the dust component of the disk. We in- 
clude the radial drift motion of solids in our model. The 
maximum radial inward drift velocity is given by 

d r P K 



2p g fi k 



(1) 



For different accretion rates, this quantity is shown in 
Fig.[2]as a function of disk location. The equations for the 
ga s pressure P<? and t he Kepler frequency £lk are given 
in IBrauer et al] ( 2008). From this maximum radial drift 
velocity, we can calculate the actual radial drift speed of a 
particle of a certain size and solve the continuity equation 
in the radial direction for all dust particle species. Radial 
diffusion due to turbulent mixing is is addition included 
in the model. 

Fig.[2]shows, that for sufficiently small accretion rates, 
grain retention will occur around the snow line close to 
3 AU inside the disk. For dust particle growth, it is even 
more important to note that there are two radial drift 
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Fig. 1. The surface gas densities for different gas accretion 
rates which we will u se in our simulations as calculated by 
Kretke &: Lira ( 20071 ) and discussed in Sec.[2j The snow line 
is located at 3 AU. 
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Fig. 2. Maximum radial drift velocity «n for different gas 
accretion rates as discussed in Sec. [2] The horizontal line 
indicates zero radial drift. Outward drift occurs for the 
two accretion rates 7 x 10 -9 and 8 x 10~ 9 M sun /yr. Hence, 
we expect dust grain retention for these two values. Note 
as well, that for these two accretion rates, there are two 
stagnation points in which the radial drift speed vanishes. 



We consider that the disk has a vertical structure. 
The vertical dust distribution is determined by turbu- 
lent diffusion and par ti cle settling towards the mid plane 
(|Dubrulle et alJ 1 199.4 ISchraoler fc Hennind l2004h . We 
solve the sedimentation/diffusion equation with a semi- 
analytical model, assuming that the particle distribution 
is in equilibrium between particle settling and turbulent 
mixing. We consider a dead zone around the midplane in 
which almost no turbulence is present. In this region, we 
adopt an a- value of 10 -5 . In the MM active regions on the 
surface of the disk , we assume an a parameter of 1.8 x 10~ 2 
( Sano et al]|l998l ). The amount of g as and dust mass tha t 
is MRI active/inactive is provided in lKretke fc Lin ( 2007 ). 
Little turbulence around t he midplane involves s mall rel- 
ative turbulent velocities (jOrmel fc Cuzzil 120071 ). Hence, 
in the radial drift stagnation points, all relative velocities 
are low and we may expect dust particles to grow to larger 
sizes and break through the fragmentation barrier in the 
midplane of the disk. 

At every point in the disk, we allow dust particles to 
coagulate and grow to larger sizes or suffer fragmentation 
due to high speed collisions. The parameter that distin- 
guishes between these two cases is the collision velocity 
of the particles. If this quantity is higher than a critical 
threshold velocity Uf, fragmentation occurs. Apart from 
the accretion rate M acc , the threshold fragmentation ve- 
locity is the second parameter that we consider. In the 
simulations, we assume five different sources of relative 
particle velocities in the disk that generate dust growth, 
namely Brownian motion, differential settling, relative 
turbulent velocities, and relative radial and azimuthal ve- 
locities. All of these sources are explained extensively in 



iBrauer et al. ( 2008 ). which also describes the numerical 
schemes adopted to solve the coagulation equation (i.e. 
Smoluchowski equation), the outcome of fragmentation, 
and the effect of cratering that is also included in the 
model. 

In the final simulations completed by IBrauer et al 



(|2008l ), particles were always smaller than the mean free 
path of the gas at any point in the disk. This is not the case 
here due to the very high gas densities. This implies that 
we must consider a different drag force regime, namely 
the Stokes regime instead of the Epstein regime. We im- 
plemented this reg ime into our model to account for the 
high gas densities (|Weidenschillindll977l ). 



stagnation points, i.e. disk locations at which the radial 
drift completely vanishes. Relative particle speeds at these 
radii are too low to induce destructive particle collisions. 
While solid particles move away from the inner diverg- 
ing point at ~ 2.9 AU, the dust tends to accumulate in 
the outer converging point at ~ 3.3 AU. It is also inter- 
esting to note here, that the azimuthal gas flow becomes 
Keplerian at these stagnation points. The differential az- 
imuthal speed between the dust and gas vanishes, which 
implies that self-induced turbulence is stablized regardless 
on how thin the dust layer becomes. 



3. Results 

Fig-Oshows the particle distribution after 1800 yrs of disk 
evolution for different accretion rates. In this simulation, 
we adopted a critical threshold fragmentation velocity of 
10 m/s. This plot indicates that particles can grow to some 
10 2 m in size around the ice evaporation front subject 
to the condition that the gas accretion rate is not too 
high. It also shows that fragmentation inhibits particle 
growth towards m-sized particles in other disk regions, 
for example around 4 AU. To unveil the importance of 
collective effects at the snow line, we estimate the dust- 
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Fig. 3. The particle distribution after 1800 yrs of disk 
evolution as discussed in Sec. [3] Shown is a contour plot of 
the surface dust density as a function of disk location and 
particle radius for four different accretion rates. The figure 
indicates that particles break through the fragmentation 
barrier if the accretion rate is not too high. The critical 
fragmentation velocity is 10 m/s in this simulation. 
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Fig. 4. As Fig. [3j but now for four different critical thesh- 
old fragmentation velocities. The figure indicates that dust 
particles can break through the fragmentation barrier if 
the theshold fragmentation velocity is at least 5 m/s. In 
this simulation, we adopted an accretion rate of M acc = 
8 x 10~ 9 M sun /yr. 



to-gas ratio in the midplane of the disk. The vertically 
integrated ratio in the case of M = 9 x 10~ 9 M sun /yr is 
of the order eo = 10/1 4 = 10~ 3 , as inferred from Figs. Q] 
and [3j The dust-to-gas ratio in the midplane is then at 
least e mid « e oy /St/a » lO^^/WJW^ = 10 > 1, which 
means that collective effects play a non-negligible role. 
Since collective effects strongly influence the radial drift 
behaviour of the dust, further investigation of this issue is 
certainly required. 

The accretion rate of a protoplanetary disk decreases 
with time. Therefore, Fig. [3] illustrates the advantage of 
our planetesimal formation mechanism at different stages 
of disk evolution. In the early stages, the gas surface den- 
sities are not significantly affected by the sharp decrease 
in dust-to-gas ratio at the snow line because the accretion 
rate is too high. Hence, we do not find that particles grow 
to become very large objects in this case. With decreas- 
ing accretion rate, the surface gas density becomes more 
and more affected by the snow line producing the forma- 
tion of large boulders as shown in Fig. [3j At this stage 
of disk evolution, planetesimal formation could occur. At 
later stages, the surface gas densities will have declined to- 
wards a level at which the entire disk becomes MRI active 
and the mechanism is unable to operate. 

The dust particle distribution after 1800 yrs, but now 
for different criticial fragmentation velocities, is shown in 
Fig.[H In this calculation, we considered an accretion rate 
of -/Vf acc = 8x 10~ 9 Msun/yr. This plot indicates that dust 
particles can overcome the fragmentation barrier and grow 
to nearly km-size if the critical threshold velocity is at least 
5 m/s. For the case Vf = 1 m/s, we do not find particles 
that are able to grow to larger than a meter in radius. 
This is due to the fact that relative turbulent velocities 



in our quiescent midplane are still of the order of several 
m/s. Even with a low residual turbulent ao-value of 10~ 6 , 
relative particle velocities would hardly decrease to below 
this value. 

The residual turbulent a -value around the midplane 
can be shifted to high er values due to t he oc currence 
of turbulent convection. iKlahr fc Henning (jl9971 ) suggest 
turbulent gas speeds of ~ 0.02c s , which correspond to ve- 
locities of ~ 16 m/s at 1 AU and an a-value of ~ 10~ 4 . 
Hence, if convection is a driving source of turbulence in 
the dis k then the resu l ts of t his Letter might change signif- 
icantly. iTurner et alj (|2007t ) showed that free charges can 
be mixed into the disk interior, producing a slight coupling 
between the midplane gas and the magnetic fields. This 
may generate an active turbulent midlane layer, which 
could trigger particle fragmentation. Further investigation 
of the influence of each of these effects is imperative. 



4. Summary and conclusions 

Even though particle fragmentation by high speed colli- 
sions is a severe obstacle to the formation of planetesi- 
mals, we have shown that dust particle growth to the size 
of enormous boulders is possible under certain circum- 
stances. If we consider dust particle coagulation in the 
presence of an evaporation front, solid particles grow to al- 
most km-sizes within only a few thousand years. The pro- 
toplanetary disk must also contain a dead zone in which 
turbulence - the main source of relative particle velocities 
- is almost entirely absent. Our conclusion is that plan- 
etesimal formation due to dust particle agglomeration is a 
possible mechanism to form large bodies in protoplanetary 
disks. 
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